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ABSTRACT 

We present 2D multi-group flux-limited diffusion magnetohydrodynamics (MHD) simulations of 
the Accretion-Induced Collapse (AIC) of a rapidly-rotating 1.92 M white dwarf, started from a 2D 
rotational equilibrium configuration. We focus on determining the dynamical role of MHD processes 
after the formation of a millisecond-period protoneutron star. We find that, compared with the 
previous study of Dessart et al. that did not include magnetic fields, magnetic stresses can lead to a 
powerful explosion with an energy of a few Bethe, rather than a weak one of at most 0.1 B, with an 
associated ejecta mass of ~0.1 M , instead of a few 0.001 M . The core is spun down by ~30% within 
500 ms after bounce, and the rotational energy extracted from the core is channeled into magnetic 
energy that generates a strong magnetically-driven wind, rather than a weak neutrino-driven wind. 
Baryon loading of the ejecta, while this wind prevails, precludes it from becoming relativistic. This 
suggests that a 7-ray burst is not expected to emerge from such AICs during the early protoneutron 
star phase, except in the unlikely event that the massive white dwarf in this AIC context has sufficient 
mass to lead to black hole formation. In addition, we predict both negligible 56 Ni-production (that 
should result in an optically-dark, adiabatically-cooled explosion) and the ejection of 0.1 M of material 
with an electron fraction of 0.1-0.2. Such pollution by neutron-rich nuclei will put strong constraints 
on the possible rate of such AICs. Moreover, being free from "fallback," such highly-magnetized 
millisecond-period protoneutron stars may, as they cool and contract, become magnet ars, and the 
magnetically-driven winds may later transition to Poynting-flux dominated and relativistic winds, 
eventually detectable as 7-Ray Bursts (GRBs) at cosmological distances. The likely low event rate of 
AICs, however, will constrain them to be, at best, a small subset of GRB and/or magnet ar progenitors. 
Subject headings: stars: white dwarfs - stars: neutron - stars: supernovae: general - neutrinos - 
rotation - Gamma-ray: bursts 



1. INTRODUCTION 

Some white dwarfs, located in binary systems, are 
thought to lead to thermonuclear runaways and Type la 
supernovae, a circumstance in which a Chandrasekhar 
mass carbon and oxygen core is incinerated, on a 
timescale of ~ls, to intermediate-mass and iron-peak 
elements. In the single-degenerate scenario, the white 
dwarf accretes material from a hydrogen or helium donor 
star. Aided by the stabilizing effect of fast and dif- 
ferential rotation (Yoon et al. 2004), it reaches the 
~1.4Mq Chandrasekhar limit and explosively burns car- 
bon and oxygen under degenerate conditions. However, 
for sufficiently high mass accretion rates, carbon and oxy- 
gen may burn under non-degenerate conditions to allow 
the formation of a high-density ONeMg core which, upon 
achieving the Chandrasekhar mass, collapses on a dy- 
namical timescale of ~100ms (Nomoto & Kondo 1991; 
Nomoto et al. 1995; Gutierrez et al. 1996; Gil-Pons & 
Garcia-Berro 2001). In the double-degenerate scenario, 
two white dwarfs in a short-period binary system eventu- 
ally coalesce to form a massive white dwarf that exceeds 
the 1^-corrected Chandrasekhar mass limit. Smooth- 
Particle-Hydrodynamics (SPH) simulations (Benz et al. 
1990; Mochkovitch & Livio 1989, 1990; Segretain et al. 
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1997; Guerrero et al. 2004; Yoon et al. 2007) predict the 
formation of a complex object composed of a cold, slowly- 
rotating core, surrounded by a hot and fast-rotating en- 
velope, and, further out, by some residual mass in a 
Keplerian disk. Depending on the temperature in the 
envelope, the mass accretion rate on the newly-formed 
white dwarf, and the mass of each white dwarf com- 
ponent, either Type la or Accretion Induced Collapse 
(AIC) are possible. Overall, the occurrence rate of the 
AIC of white dwarfs is difficult to determine reliably, al- 
though it seems unlikely that they occur more frequently 
than once per 20-50 Type la event (Yungelson & Livio 
1998,2000; Madau et al. 1998; Fryer et al. 1999; Nele- 
mans et al. 2001ab; Blanc et al. 2004; Mannucci et al. 
2005; Belczynski et al. 2005; Greggio 2005; Scannapieco 
& Bildsten 2005; Dessart et al. 2006a, hereafter D06). 

In this work, we focus on the subset of accreting white 
dwarfs that form ONeMg cores, which, due to their high 
central density and mass, collapse to form a neutron star. 
The formation path followed, either through accretion 
from, or coalescence with, a companion, suggests a large 
angular momentum budget for the resulting white dwarf. 
Such white dwarfs rotate differentially, with a specific an- 
gular momentum that increases outward, and can be dy- 
namically stable well above the canonical Chandrasekhar 
mass of ~1.4Mq, reaching theoretical masses as high 
as 4.1 M (Ostriker & Bodenheimer 1968; Hachisu 1986; 
Yoon & Langer 2005). The precise distribution of the an- 
gular momentum inside the star is not accurately known, 
in particular since magnetic torques may spin down the 
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core and enforce solid-body rotation, producing an even 
more slowly rotating core than presently predicted by 
current SPH simulations (~20s period). In any case, 
the overlying envelope is fast-rotating and the resulting 
progenitor is strongly aspherical (Guerrero et al. 2004; 
Yoon et al. 2007) 3 . 

The collapse of ONeMg cores has been investigated in 
the past both in the context of white dwarfs (Baron et 
al. 1987b; Woosley & Baron 1992; Fryer et al. 1999; 
D06) and in the context of 8-10 M stars (i.e., low-mass 
massive stars; Hillebrandt et al. 1984; Mayle & Wilson 
1988; Baron et al. 1987a; Kitaura et al. 2006; Bur- 
rows et al. 2007c). In D06, we presented results of ra- 
diation hydrodynamics simulations of the AIC of white 
dwarfs that extended these past investigations simul- 
taneously to include multi-dimensionality, multi-group 
flux-limited diffusion for multi- flavor neutrino transport, 
and, perhaps most importantly, rapid rotation. We fo- 
cused on the 1.46 M and the 1.92 M models of Yoon 
& Langer (2005) and found for the latter that fast ro- 
tation led to numerous striking features, e.g., aspher- 
ical collapse, aspherical bounce, polar-confined explo- 
sions, polar-enhanced neutrino radiation fields, gravity- 
darkened neutron stars, oblateness of the neutron stars 
formed, and a residual Keplerian disk. In this, the dif- 
ferential rotation that results from core collapse sets the 
stage whereby, from initial seeds in the poloidal com- 
ponent, a magnetic field may exponentially grow on 
a timescale of a rotation period by the magnetorota- 
tional instability (MRI; Balbus & Hawley 1991; Pessah 
& Psaltis 2005,2006ab; Etienne et al. 2006), acting on 
the poloidal and on the toroidal fields. If equipartition 
between the magnetic energy and the differential rota- 
tional energy is reached at saturation, very substantial 
fields may be generated for the fastest rotators (Shibata 
et al. 2006). The 1.92 M AIC white dwarf studied in 
D06, with an initial T/\W\ of 10% (W is the gravitational 
energy of the initial equilibrium configuration, and T is 
the initial rotational kinetic energy), represents a good 
candidate for MHD effects after core bounce. Unfortu- 
nately, it is not feasible, even with modern supercom- 
puters, both to follow the multi-dimensional evolution 
of a multi-group neutrino field for three species and to 
resolve finely enough the scales over which the magne- 
torotational instability operates. To circumvent this dif- 
ficulty, we start with fields that are strong enough so 
as to reach soon after core bounce (after a few tens of 
rotation periods) the saturation values they would have 
reached had the MRI operated. Magnetic- field ampli- 
fication then results exclusively from compression and 
winding. This is the approach followed in Burrows et 

3 Observationally-inferred rotation rates of white dwarfs are, 
however, generally modest, a fact attributed to the copious an- 
gular momentum losses during their evolution, even for isolated 
white dwarfs (Koester et al. 1998; Berger et al. 2005; Valyavin et 
al. 2005; Karl et al. 2005). Magnetic, fast-rotating white dwarfs in 
binary systems may also undergo considerable spin-down, caused 
by magnetic torquing of the envelope (Ikhsanov 1999). Merger 
events show, by contrast, a potential for the production of fast- 
rotating white dwarfs, with properties more germane to the present 
discussion (Ferrario & Wickramasinghe 2005). If fast rotation of 
either the core or the envelope were compromised, the outcome 
of the AIC would be more suitably described by the simulation of 
D06 for the 1.46 M model, rather than that for the 1.92 M model 
described here or in D06. 



al. (2007b), where we presented results for magnetically- 
driven supernova explosions of massive star progenitors 
in the context of fast rotation. 

The central motive for this work is to go beyond the 
investigations of D06 and to study the potential role of 
MHD effects, adopting initial fields that will be commen- 
surate at saturation with values expected from equiparti- 
tion with the free energy of rotation in the differentially- 
rotating surface layers of the rotating collapsed core. 
Thus, we investigate the MHD effects that may arise dur- 
ing the collapse phase and the few hundreds of millisec- 
onds that directly follow the formation of the neutron 
star. Guided by our simulations, we comment on past 
and current suggestions that the AIC of white dwarfs 
are potential Gamma- Ray Burst (GRB) and/or magne- 
tar progenitors, or are important production sites of r- 
process nuclei. 

In the next section, we discuss (briefly) the VUL- 
CAN/2D code, and the numerical setup, as well as the 
1.92 M AIC model we focus on. The details of the code 
can be found in Livne et al. (2004,2007), Dessart et 
al. (2006b), and Burrows et al. (2007a), and the AIC 
initial model properties are given in D06. In §3, we de- 
scribe the results of our simulations, and compare a ref- 
erence model with weak initial fields to another model 
with strong initial fields. We discuss our results in the 
broader astrophysical context in §4. In §4.1, we sum- 
marize the main MHD effects we identify, and in §4.2 
we discuss the potential nucleosynthetic yields. In §4.3, 
we discuss potential electromagnetic displays from AICs, 
and whether they can be 7-ray bursts, and in §4.4 we dis- 
cuss the ultimate fate of these rare objects. Finally, in 
§5, we present our conclusions. 

2. MODEL AND PROGENITOR PROPERTIES 

The models presented in this work derive from the 
1.92 M model of D06 (see their §2), originally produced 
in Yoon & Langer (2005). The initial density, tempera- 
ture, composition, and angular velocity distributions are 
coincident with those used in D06. To summarize, the 
central density is p c — 5 x 10 10 gcm -3 , the temperature 
is given by T(r,z) = T c (p c /p(r, z)) - 35 with T c = 1.3 
x 10 10 K (r and z are the standard cylindrical coordi- 
nates), the initial electron fraction (Y e ) is 0.5 throughout 
the grid, and the central angular velocity is 20.4 rads -1 
(with a period of 0.3 s). To minimize the dynamical influ- 
ence of the medium surrounding the white dwarf, which 
needs to be initially defined in our Eulerian approach, we 
assign it properties that correspond to the low-density 
low-temperature edge of our equation of state (EOS, see 
also below) table, i.e., a density of 10 3 gem -3 , and a tem- 
perature of 4 x 10 8 K. Furthermore, we also assign it a 
negligible angular velocity. Initially, the polar (equato- 
rial) radius is 660 km (2350 km), the total angular mo- 
mentum is 0.1092 B • s (1 B=10 51 erg), the total rotational 
energy T is 1.057 B, and the gravitational energy \W\ is 
12.69 B. This model thus starts with a ratio T/\W\ of 
0.083, indicative of the large angular momentum budget 
in the progenitor white dwarf. We refer the reader to §2 
and Fig. 1 of D06 for a detailed description of this initial 
model, with 2D representations of the original distribu- 
tion of the density, temperature, and angular velocity. 

The magnetic field morphology we adopt is uniform 
inside a sphere of radius 600 km (the poloidal compo- 
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nent has a z-orientation, i.e., is aligned with the axis 
of rotation), and dipolar beyond. Note that the origi- 
nal model of Yoon & Langer (2005) corresponds to an 
equilibrium configuration that results from the balance 
between gravitational, rotational, and thermal energy, 
and conveys no information on the magnetic-field mor- 
phology and magnitude. One may argue that magnetic 
torques would enforce solid-body rotation in the initial 
model we adopt from Yoon & Langer (2005), making our 
approach inadequate. We propose in this work that even 
such super- Chandrasekhar white dwarfs are unlikely to 
lead to black hole formation, so the differential rotation 
of our initial model does not compromise our conclusions 
in this respect. Our main assumption is that such torques 
may alter the rotational profile, but not reduce signifi- 
cantly the overall rotational energy budget of the white 
dwarf, preserving its status as a fast rotator. Further- 
more, the initial differential nature of the rotation is not 
essential since gravitational collapse will naturally lead 
to it. For the binary merger scenario, magnetic torques 
would likely not have sufficient time to sizeably torque 
the white dwarf before core collapse. Little is known 
about the core magnetic fields of white dwarfs, and even 
less on that fraction of them that may lead to AIC and ro- 
tate fast, nor do we know how the poloidal and toroidal 
components relate. Again, the key here is that we as- 
sume that starting from a given initial seed poloidal field, 
the magnetorotational instability will result in the expo- 
nential growth of the magnetic energy after bounce, and 
then saturate, after a few tens of rotation periods, at 
the field we obtain here. By that time, we expect that 
the system will have lost memory of the magnetic field 
distribution it started with. In practice, we perform sim- 
ulations with an initial poloidal field of 10 11 G for model 
Bll and of 10 12 G for model B12, while in both cases we 
use a toroidal field of 10 8 G (these values are at the white 
dwarf center). Note that the magnetic field amplification 
at bounce resulting from compression alone will lead to 
an enhancement of the magnetic field by a factor of only 
(Pnuciear/Pc)^ 3 ^ °° 300, because of the very high initial 
central density (p c = 5 x 10 10 gem -3 ). For a star with an 
initial central density of 5 x 10 9 g cm -3 , a typical value at 
the onset of core collapse for a massive-star progenitor, 
the field amplification would be by a factor of ~1200 in- 
stead, so that an initial field of 2 x 10 11 G would achieve 
a comparable magnetic field at bounce. Hence, initial 
field values must always be interpreted together with the 
initial central densities to properly gauging the amplifica- 
tion generated by compression. The choice of a predom- 
inantly poloidal component is essential in our approach, 
allowing for the growth of the toroidal component after 
the initial compression phase. If the MRI operated, the 
exponential growth of both components of the magnetic 
field, permitted by the large budget of free energy of ro- 
tation, would reach saturation values that are orders of 
magnitude above the initial seed values. Hence, in this 
context, the main limiting factor is the rotational energy 
budget rather than the initial properties of the field. The 
evolution of the Bll model is similar to that of the non- 
magnetic 1.92 M AIC model presented in D06, and we 
use it to gauge the importance of MHD effects in the B12 
model. The lack of magnetic effecs in the Bll model is 
a clear sign that we do not resolve the MRI, since oth- 



erwise we would observe magnetic fields that saturate at 
values similar to those seen in the B12 model, and we do 
not. Note that the 1.92 M progenitor model we discuss 
in this work is not meant to be representative of AICs in 
Nature, as, indeed, these may come with a range of an- 
gular momenta or masses. But our B12 model is chosen 
because it is one example where all conditions are prime 
for the generation of MHD effects, thus occupying one ex- 
treme when compared to, for example, the more slowly 
rotating and quasi-spherical 1.46 M model simulated in 
D06. This approach allows us, we hope, to bracket the 
range of outcomes for AICs. 

We adopt the same spatial grid as D06, with a Carte- 
sian mapping of the region inside a spherical radius of 
20 km and a spherical-polar grid out to a maximum ra- 
dius of 4000 km. The maximum resolution of 0.48 km is 
reached in the Cartesian inner region, and the minimum 
resolution is 100 km at the maximum radius, with 71 
regularly-spaced angular zones covering 90°. In D06, we 
found that the top/bottom symmetry was perfectly pre- 
served in the 180° simulations of the 1.46 M AIC model. 
We speculate that the post-bounce time to explosion is so 
short that convection and the Standing- Accretion- Shock- 
Instability (SASI; Blondin et al. 2003; Foglizzo et al. 
2006) cannot grow to any sizable level, perhaps inhib- 
ited by the strong differential rotation (Burrows et al. 
2007b). Hence, we save time and opt for higher resolu- 
tion by restricting our investigation to one hemisphere. 

The simulations we present here constitute the state- 
of-the-art, combining in the dynamical equations both 
the effects of magnetic field and neutrino transport. As 
in D06, we employ the approximation of Mult i- Group 
Flux-Limited Diffusion, using 16 logarithmically-spaced 
energy groups between 2.5 MeV and 220 MeV, for three 
neutrino species (electron neutrino z/ e , anti-electron neu- 
trino P e , and a third entity grouping the remaining neu- 
trinos, designated We also employ a tabulated 
version of the Shen (1998) EOS (see Dessart et al. 2006b 
for details). Finally, the gravitational potential is com- 
puted using a multipole solver, with a decomposition 
in spherical harmonics that extends to I = 32. Again, 
the VULCAN/2D code and our standard approach to 
modeling core collapse are described in Livne et al. 
(2004,2007), Dessart et al. (2006b), and Burrows et al. 
(2007a), while further details on the AIC of white dwarfs 
have been given in D06. 

3. DESCRIPTION OF THE EXPLOSION 

In this section, we present the collapse, the bounce, 
and the explosion phases of the 1.92 M AIC white dwarf 
model computed by, and with the method of, Yoon & 
Langer (2005). In contrast to a similar investigation per- 
formed and analyzed in great detail in D06, magnetic 
fields are now included, leading to little qualitative, but 
considerable quantitative, changes in the results. For 
completeness, we briefly review the important phases of 
the explosion arising from such AIC of white dwarfs. Im- 
portant results are discussed again, and in a broader con- 
text, in JH 

Due to the high central density of such white dwarfs 
(above ~10 10 gcm -3 , and here set to 5 x 10 10 gcm -3 ), 
electron captures on nuclei lead to the collapse of the 
ONeMg core on a near dynamical timescale. About 37 ms 
after the start of the simulation, the inner core reaches 
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Fig. 1. — Top: Colormap of the entropy for the Bll (left) and B12 (right) models, at 816 ms and 670 ms after bounce, respectively, and 
in the inner 1000 x 1000 km 2 , using black vectors to represent the velocity field (the length is saturated at 5000 kms -1 , corresponding to 
5% of the width of the display). We also overplot isodensity contours, shown for every decade starting at 10 14 gcm -3 . Bottom: Same as 
for the top panels, but extending the maximum displayed radius to 4000 km, the maximum grid radius in our simulations. The saturation 
length for the velocity vectors is now 10000 kms -1 . (See text for discussion.) 



nuclear densities, i.e., >2 x 10 14 gcm -3 . A shock wave 
is born that propagates outward, but weakens due to 
the huge energy losses associated with the burst of elec- 
tron neutrinos and the photodissociation of nuclei in the 
infalling outer-core material. Along the equator, where 
the accreting material enjoys sizable centrifugal support, 
the shock first propagates fast due to the reduced ac- 
cretion rate, but eventually stalls due to its sustained 
magnitude. By contrast, along the polar direction where 
the density profile is much steeper, after at first stalling 
due to the higher accretion rate, the shock soon encoun- 
ters the layers that were formerly at the surface of the 
white dwarf. As no further accretion occurs, the shock 
blasts out of the white dwarf along the pole only 80- 
100 ms after bounce and an explosion is initiated. Up to 
that point, the simulations Bll and B12 show the same 
qualitative and quantitative properties, suggesting that 
magnetic effects on the dynamics are weak or negligible 



during the first ~100 ms post-bounce evolutionary phases 
of the AIC. 

However, it takes ~100ms in the B12 model to amplify 
the magnetic field to a value such that magnetic pressure 
and gas pressure are comparable, at first in the polar re- 
gions at a few tens of kilometers from the protoneutron 
star (PNS) center. Compression from the progenitor to a 
high-density collapsed configuration is the early cause of 
magnetic field amplification, but later it arises from the 
winding of the poloidal field and the continued accretion 
of the magnetized envelope, spinning up as it collapses 
(see discussion in Burrows et al. 2007b). The toroidal 
field eventually dominates over the poloidal field compo- 
nent in the inner few hundred kilometers by, typically, 1- 
2 orders of magnitude. We speculate that the magnetic 
properties at ~100ms after bounce are comparable to 
what would obtain with a fully resolved MRI. Note that 
in our simulations, no significant convection develops be- 
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Fig. 2. — Left: Time evolution for the 1.92 Mq model of Yoon & Langer (2005), with either low (black) or high (red) initial magnetic 
fields, of the total energy associated with the ejecta (solid line; accounting for gravitational, radial and rotational kinetic, thermal, and 
magnetic energy components) and the total integrated net neutrino energy gain (broken line). Right: Same as left, but for the Bernoulli 
(solid line) and Poynting (broken line) powers in the ejecta, computed by integrating the corresponding flux over a radial shell at 500 km. 



hind the shock, as required for the a 2 and the a-fl dy- 
namo invoked for magnetic-field amplification by Duncan 
& Thompson (1992) and Thompson & Duncan (1995). In 
any case, in the B12 model after ~100ms, the resulting 
large magnetic pressure and associated gradient at the 
surface of the neutron star power a strong magnetically- 
driven outflow that expands in the previously excavated 
polar regions of the progenitor white dwarf. There is here 
an analogy with magnetic tower formation in magnetic, 
differentially-rotating disks (Lynden-Bell 1996, 2003; Uz- 
denski & MacFadyen 2006), since the general picture is 
one of magnetic field amplification in a medium that ex- 
erts pressure anisotropically. Here, the presence of exca- 
vated polar regions and a dense centrifugally-supported 
fat disk at lower latitudes offers a natural confinement 
mechanism for the magnetically-driven explosion. Hoop 
stresses act poorly to confine the blast, since wind parti- 
cle trajectories originate at the neutron star surface over 
all latitudes between ~40° and 90° 4 . By comparison, 
in the Bll model, there is a "quiescent" stage after the 
original blast, with some fallback of a fraction of the 
ejecta. The explosion is revived after ~300ms by the 
birth of a neutrino-driven wind, as discussed in detail 
in D06. This wind, and the neutrino energy deposition 
that drives it, are the source of the energy of the explo- 
sion in the context of low or negligible magnetic fields. 
At high B-fields, the neutrinos play only a secondary, 
sub-dominant, role, and then only after a few hundred 
milliseconds after bounce. 

We have carried out the Bll and B12 simulations up to 
816 ms and 670 ms, respectively, after bounce. This con- 
stitutes only the very early evolution of the cooling PNS, 
which typically takes tens of seconds in total when in iso- 
lation (Burrows & Lattimer 1986). We show the entropy 
(colormaps) and density (contour lines) distributions in 
Fig. [T] for the Bll (left column) and B12 (right column) 
models at the end of each simulation, encompassing the 
full grid (bottom row) or the inner region (top row). We 
can identify three different regions: an inner region where 
an oblate neutron star resides, a prolate mass distribu- 
tion in the ejecta that is mostly confined to polar or 

4 The analogy to an infinitely long thin disk threaded by a mag- 
netic field is not strictly applicable. 



high-latitude regions (above ~40°), and a fat and Keple- 
rian disk that closely resembles the original lobes of the 
progenitor white dwarf. We now review in detail each 
component, together with the global energetics of each 
explosion. 

In the left panel of Fig. [2j we show the time evo- 
lution, for both models (B12: Red; Bll: Black), of 
the explosion energy (solid line) and integrated net neu- 
trino energy gain as a function of time after bounce (we 
limit this integration to regions with a mass density less 
than 10 10 gem -3 ; broken line). With high initial B-field, 
the explosion energy reaches ~1 B, two orders of magni- 
tude higher than in the corresponding low B-field case 
(which explodes under-energetically, but easily, nonethe- 
less), and one order of magnitude larger than the inte- 
grated net gain through neutrino energy deposition. At 
low B-field, the neutrinos play a critical role in power- 
ing the explosion, but the explosion is symptomatically 
weak, on the order of 0.01 B. In the right panel of Fig|2j 
we show the power injected into the ejecta, associated 
both with the Bernoulli 5 flux (solid line) and the Poynt- 
ing flux (broken line) integrated through a sphere at 
500 km. Although the Poynting flux decreases roughly 
as the inverse square of the radius above the neutron 
star surface, it is larger in the B12 simulation than the 
Bernoulli flux in the Bll simulation. Combined with 
the strong Bernoulli flux, the explosion power reaches 
-2BS- 1 . 

This gain in explosion energy stems partly from the 
extraction of the free energy of rotation of the PNS, the 
difference between the kinetic energy of rotation and the 
corresponding energy for solid-body rotation at the same 
total PNS angular momentum. Selecting the regions in- 
side of a density cut of 10 8 , 10 10 , or 10 12 gem -3 , the total 
rotational energy difference between the Bll and the B12 
cases is on the order of ~10B. While a large fraction of 
this energy is used to expand the PNS and, therefore, 

5 The Bernoulli flux is defined by the quantity pVr x (e t h + 
e kin + e rot + P/p — GMpns/ R), where Vr is the radial velocity, 
eth is the thermal energy, e^m and e ro t are the (r, z) and rotational 
kinetic energies, G the gravitational constant, Mp^s the PNS mass 
(taken as the cumulative mass of material with a density greater 
than 10 10 gem -3 ), and R the spherical radius. Both the Bernoulli 
flux and the Poynting flux are then integrated through a sphere at 
500 km to give a power, or a luminosity (Fig. [2)- 
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Fig. 3. — Top: Colormap of the logarithm of the ratio of magnetic pressure to gas pressure for the Bll model (left) at 816 ms after 
bounce, and for the B12 model (right) at 670 ms after bounce, using black vectors to represent the Poynting luminosity (the Poynting flux 
is multiplied by a factor 4ttK 2 and the resulting length is saturated at 0.1 Bs _1 , which corresponds to 5% of the width of the display). We 



the neutron star vicinity to prolate at a few thousand kilometers, 
kilometers. The isodensity contours are now plotted up to 10 14 gcm 

do work against gravity (where it is the most costly, i.e., 
at the base of the potential well), some of this energy 
is made available to the explosion and the growth of the 
magnetic field, whose associated gradient at the PNS sur- 
face is responsible for driving mass outward at the escape 
speed (Shibata et al. 2006; Burrows et al. 2007b). In 
the B12 simulation, the PNS density structure is signifi- 
cantly modified, with isodensity contours more spatially 
separated, and with less contrast between the polar and 
equatorial radii (see the top row of Fig. [1] for an illus- 
tration). Magnetic support is partly responsible for the 
more extended PNS configuration. As a result, the PNS 
in the B12 model suffers significant spin-down, another 
indication that rotational energy has been extracted from 
the compact object by magnetic torques. The average 
period in the inner 10 10 gcm -3 increases from 15 ms to 
21ms, and in the inner 10 12 gcm -3 increases from 2 ms 



Bottom: Same as top, but now zooming in on the inner few hundred 
3 . (See text for discussion.) 

to 2.5 ms, for models Bll and B12, respectively. 

In Fig. El we illustrate the magnetic properties at the 
last computed time in each model, for both large and 
small scales, by showing colormaps of the ratio of the 
magnetic pressure to the gas pressure. These are over- 
plotted with vectors depicting the Poynting luminosity 
(flux scaled by the quantity 4ttR 2 , where R is the spher- 
ical radius). Note how in the B12 model the magnetic 
pressure is large and comparable to the gas pressure in 
most of the ejecta. The ejecta retain a higher Poynting 
flux along the polar direction, even at large distances of a 
few thousand kilometers, and this is likely caused by the 
significant confinement of the inner jet. At mid-latitudes, 
the Poynting flux drops faster than the inverse square of 
the radius. The only excluded regions, where the mag- 
netic field has no substantial effect, are the neutron star 
core (layers with a local density greater than 10 12 g cm -3 ) 
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Fig. 4. — Left: Colormap of the logarithm of the poloidal field magnitude for the B12 model at 670 ms after bounce, using black vectors 
to represent the velocity (their length is saturated to 7% of the width of the diplay for a magnitude of 10000 kms ). We also overplot 
isodensity contours, for every decade starting at 10 14 gcm -3 Right: Same as left, but for the logarithm of the toroidal field magnitude. 
(See text for discussion.) 

and the lobes of the progenitor white dwarf (low-density, 
non-collapsed regions of the progenitor white dwarf). In 
these latter regions, characterized by large vorticity but 
little bulk radial advection, we notice large fluctuations 
in the magnetic to gas pressure ratio, stemming mostly 
from fluctuations in the magnetic field energy. Although 
we cannot identify strong convection, stabilized by the 
differential rotation, moderate velocity shear can induce 
such local enhancements of the magnetic field. In Fig. [H 
we present the spatial distribution in the inner 500 km 
of each of the poloidal (left) and toroidal (right) compo- 
nents of the magnetic field for the B12 model at 670 ms 
after bounce. Unable to resolve the MRI, we predict the 
dominance of the toroidal component of the magnetic 
field over the poloidal component, typically by about two 
orders of magnitude, with the exception of the axis re- 
gion where it is the strong radial advection of material 
that generates a strong poloidal field from the original 
strong toroidal field of the corresponding mass parcel 
(in the Bll model, equipartition is reached in this re- 
gion, but the gas density and pressure are two orders 
of magnitude lower than in the same region in the B12 
model). While the relative poloidal and toroidal field 
magnitudes produced by an operating MRI would likely 
differ, the total magnetic energy we obtain may agree 
approximately, as it is fundamentally extracted from the 
(given) free energy of rotation available in the differen- 
tially rotating surface layers of the PNS. As in Fig.[3l we 
show in Fig. [5] the ratio of the magnetic pressure to the 
gas pressure, but at 10 ms, 50 ms, 100 ms, and 670 ms 
after bounce. This figure portrays how the pressure is 
more and more magnetically-dominated directly outside 
the neutron star; the magnetically-driven wind advects 
highly- magnetized material to large distances. 

Mass loss in the explosion is considerably enhanced in 
the B12 model, the cumulative mass expelled being two 
orders of magnitude larger than in the Bll model, i.e., 
0.1 M compared with a few x 0.001 M (Fig. [6]). In both 
cases, accretion in the equatorial regions, where material 



is centrifugally supported, is of the order of 0.01 M s -1 . 
In the low B-field case, this is nearly compensated by 
the mass lost through the neutrino-driven wind. In the 
B12 model, the mass loss rate through the magnetically- 
driven wind is nearly two orders of magnitude larger, at 
0.7 M s -1 , and the nascent neutron star actually loses 
mass. Interior to a density cut of 10 10 gcm -3 , the PNS 
baryonic mass is 1.32 M for the B12 model, compared 
with 1.47 M for the Bll model (both at -500ms af- 
ter bounce). At the end of the B12 simulation and at 
the maximum radius of 4000 km, the ejecta density is 
10 5 -10 6 gcm -3 and its velocity is 40000 kms" 1 . This is 
about two orders of magnitude denser and 20% slower 
than the neutrino-driven wind we obtain at comparable 
times in the Bll model. Hence, while the ejecta are mag- 
netically driven, they are not Poynting-flux dominated 
and they are even more non-relativistic. 

In D06, we found that the electron fraction of the 
ejected material was bimodal, with symmetric material 
of high entropy (i.e., — 25 &b baryon -1 ) ejected along the 
poles, and neutron-rich material of lower entropy (i.e., 
—15 &b baryon -1 ) at mid-latitudes. Here, the electron 
fraction is clustered around a neutron-rich mean, with 
Y e = 0.25 and an entropy of —8 &b baryon -1 . We il- 
lustrate these results in Fig. [7l with histograms for the 
electron fraction and entropy distributions of the ejecta 
material. Given the — 0.1 M of material expelled in the 
B12 model with such characteristics, these yields may 
have some significance for the composition of the inter- 
stellar medium in r-process nuclei (see §4.2p . as suggested 
by Wheeler et al. (1998). 

The long-term development of the explosion is influ- 
enced by the early excavation of the polar regions of the 
white dwarf. Hence, the explosion remains confined to 
the poles, extending only down to a —40° latitude. Some 
fraction of the material on the polar-facing side of the 
excavated white dwarf is entrained by the ejecta, the ve- 
locity shear at the interface inducing Kelvin-Helmholtz 
rolls that mix low and high Y e material (the lower reso- 
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Fig. 5. — Colormap of the logarithm of the ratio of magnetic pressure to gas pressure for the B12 model, at 10 ms (top left), 50 ms (top 
right), 100ms (bottom left), and 670ms (bottom right) after bounce, using white isodensity contours starting at 10 14 gcm -3 shown for 
every decade in density, as well as black vectors to represent the velocity, saturated at a length of 7% of the width of the display for a 
magnitude of 5000 km s -1 . 

lution of our simulations in the region of interest, around 
a few hundred kilometers, underresolves this instability, 
but some rolls are clearly seen). At the end of the sim- 
ulation, this interface has visibly been eroded and the 
lobes of the white dwarf progenitor subtend a smaller 
opening angle (compare the left and right panels in the 
top row of Fig. [1]). Interestingly, most of this entrained 
material does not reach the escape speed from this AIC 
progenitor, of the order of 15000 kms -1 , and ultimately 
remains bound to the system. In practice, its trajectory 
curves back towards the equatorial plane, which it hits at 
a radius of 2000-3000 km from the PNS. Hence, material 
with low Y e is being fed into the backside of the disk, 
from regions that were originally near the pole 6 . If the 

6 Note that, initially, in our Eulerian approach, we fill the am- 
bient medium around the progenitor white dwarf with low-density, 
low-temperature material, which, nonetheless, may do work against 
the eject a expansion and lead to the fall back of some material at 
the outskirts of the disk. At the start of the simulation, this ambi- 
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Fig. 6. Time evolution for the Bll (black) and B12 (red) 
models of the instantaneous integrated mass flux accreting (solid) 
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Fig. 7. — Top: Colormap of the electron fraction Y e at the end of the simulation for model Bll (left) and model B12 (right), with white 
isodensity contours for the density shown every decade starting at 10 10 gem -3 , and white velocity contours saturated at 7% of the width 
of the display for a magnitude of 10000 km s _1 . Bottom: Three-dimensional histogram of the distribution of the ejecta mass as a function 
of electron fraction, Y e , and entropy, S, for the Bll (left) and B12 (right) models. Notice the larger total ejecta mass, and the stretching 
of the electron-fraction distribution to values as low as 0.1, in the B12 model. (See text for discussion.) 



explosion persists over a long time, such ablation may 
completely erode and destroy the original lobes of the 
progenitor white dwarf, although perhaps it may merely 
redistribute this material to larger distances still within 
the white dwarf potential well. Material located beyond 
a distance of ~1000km from the PNS and above a ~40° 
latitude follows an essentially radial, ballistic, trajectory 
and escapes the gravitational potential of the massive 
white dwarf. 

At the end of both simulations, the neutron stars 
have a strongly oblate structure and rotate fast (see 
also the comparable equilibrium configurations of Liu & 
Lindblom 2001). Although the pole-to-equator radius 
contrast is essentially preserved, in the B12 model the 

ent pressure P a mbient is about 8x 10 19 dyne cm -2 , over four order 
of magnitude smaller than the pressure in the surface layers of 
the white dwarf. The PdV work done by the ejecta to expel this 
ambient material translates into an energy penalty on the order 

Of ^ambient X ( |^^max 7T#eq X 2 ^pol) ~ 2 X 10 46 erg, where 

i^max is the maximum radius of the grid, and .Req and .Rpoi are 
the equatorial and polar radii of the progenitor white dwarf. This 
energy penalty is small in comparison with the final energy of the 
explosion. 



10 10 gcm -3 isodensity contour, for example, is stretched 
out to ~100km at mid-latitudes, and to only ~50km in 
the Bll model (compare the contour morphology in the 
left and right panels, bottom row, of Fig. [I]). The same 
holds for the temperature, which remains large along the 
pole even at 100 km (~3MeV in the B12 model com- 
pared with ~1 MeV in the Bll model). Higher magnetic 
fields modify the hydrostatic configuration of the neu- 
tron star, leading to its spatial extension, and the subse- 
quent peeling back of the corresponding regions by the 
magnetically-driven wind. 

The post-bounce neutrino signals of the AIC of the 
1.92 M white dwarf we have modeled are strongly 
latitude-dependent 7 with a systematic increase (irrespec- 
tive of neutrino type) by up to a factor of two from the 
equatorial to the polar directions (D06). The gravity- 
darkened fast-rotating PNS has a much lower tempera- 
ture along the equator, which inhibits the emission, es- 
sentially thermal, of V M " neutrinos, and the flatter den- 

7 Note that there was little neutrino field anisotropy for the 
1.46 Mq model described in D06, due to the more spherical shape 
of the newly- formed neutron star. 
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Fig. 8. — Left: Time evolution of the electron neutrino (solid line), the anti-electron neutrino (dashed line), and the "z^" neutrino 
luminosity (dash-dotted line) along the equatorial (black) and polar (red) directions, for model Bll, at R =400 km, and starting 100 ms 
after bounce. Here, neutrino luminosities are integrated over neutrino energy, and correspond to the quantity 4tt R? F v (t , R, 0) , adopting 
the flux along the polar direction (0 = 0) or along the equatorial direction (0 = 90°). Right: Same as left, but for model B12. Notice the 
systematic reduction of all neutrino fluxes along the equatorial direction (neutrinos fluxes decrease in a continuous fashion towards lower 
latitudes), an effect that is essentially geometrical and stems from the larger angle subtended by the PNS as seen from higher latitudes. 
This effect is primarily caused by the fast rotation of the PNS, rather than its magnetic properties. Notice also that, for a given time after 
bounce, the electron neutrino luminosity is reduced in the B12 model, and that the anti-electron neutrino and the "^"neutrino luminosities 
remain comparable to their Bll-model counterparts. (See text for discussion.) 



sity profile along the equator modifies capture rates. The 
enhancement in neutrino luminosity at higher latitudes, 
whatever the species, is a purely geometrical effect and 
results from the increased subtended angle occupied by 
the PNS for viewing angles at higher latitudes (see also 
Walder et al. 2005). Irrespective of the energy group and 
the neutrino species, the corresponding "radiating" neu- 
trinosphere is oblate and occupies the largest solid-angle 
for observers along the poles (see D06 for further discus- 
sion on this latitudinal variation of neutrino fluxes). In 
addition, we note that, for a given time after bounce, 
the electron neutrino luminosity is reduced in the B12 
model compared to the Bll model, by ~30% along the 
equator and ~40% along the pole, while the luminosi- 
ties associated with the u e and the neutrinos match 
their counterparts in the two models within a few per- 
cent. Note that, provided it is not later ejected, the de- 
layed mass accretion, over viscous timescales (see 
of the lobes of the progenitor white dwarf has the poten- 
tial to power the neutrino luminosity at a low level over 
a longer time than would typically obtain for isolated 
neutron stars. 

4. DISCUSSION 

4.1. Magnetic effects 

The previous section highlighted the main features and 
phases of the AIC of white dwarfs, for both weak and 
strong magnetic fields. The main effect of magnetic fields 
is to trigger a much more powerful explosion, here reach- 
ing ~1B after just 300 ms after bounce, tapping some 
free energy of rotation of the PNS, and ejecting mass 
at a rate of a few x 0.1 M s _1 . We thus find that if 
magnetic pressure becomes comparable to gas pressure 
at the neutron star surface, a magnetically-driven wind 
is initiated whose properties are boosted when compared 
to the pure neutrino-driven wind case (i.e., with weak or 
negligible magnetic field), the latter with mass loss rates 
of a few x 0.001 M s _1 and global ejecta energies of at 
most ^0.1 B. Note, in particular, that the free energy of 
rotation is also available at the same level in the weak 
magnetic- field case (model Bll); it is only that the means 



to extract it are insufficient. Magnetically- aided, the ex- 
plosion is sustained at a much higher level, rendering the 
neutrino component sub-dominant (Fig. [2]). 

The magnetic contribution to the pressure balance at a 
few tens of kilometers can be significant and can deform 
the PNS compared to the weak field case, sphericizing 
the isodensity contours (see, e.g., the bottom panels of 
Fig.[T]). The energy to amplify the magnetic field, to puff- 
up the neutron star (although mostly along the poles), 
and to drive a more powerful explosion stems from the 
tapping of core rotational energy (where it is the most 
abundant), which is also manifest in the spin down of the 
core, in analogy with the extraction of angular momen- 
tum in accretion disks (Blandford & Payne 1982). By 
~500ms after bounce, the average PNS period 8 is 30% 
larger in the high magnetic field case, and is ~20ms. 

4.2. Nucleo synthetic yields 

Guided by the results of Hillebrandt et al. (1984), 
Wheeler et al. (1998) suggested promptly exploding 
ONeMg cores as the production sites of r-process nu- 
clei, given the low electron fraction (Y e ~0.2) and the 
low entropy of their ejecta (S ~ 15 &b baryon -1 ; see also 
Qian & Woosley 1996; Hoffman et al. 1997). However, 
radiation-hydrodynamics simulations of the collapse of 
Fe or ONeMg cores demonstrated that the explosion can 
never be prompt (can never be on timescales of a few 
milliseconds) and that, instead, the shock stalls. Even- 
tually, in the context of the ONeMg cores of low-mass 
massive stars or AICs, the stalled shock overcomes ac- 
cretion, because of the extremely steep density gradient 
above the core (Kitaura et al. 2006; D06; Burrows et al. 
2007ac). 

The results of our present simulations confirm that 
there is no prompt explosion in the AIC of white dwarfs 
and that, indeed, ejecta with low electron fraction are 
produced. At low magnetic field, the ejecta are primar- 
ily at Y e values between 0.25 and 0.5, but with a total 

8 The average PNS period is computed assuming solid-body ro- 
tation for the same angular momentum budget in the PNS, with 
densities greater than 10 10 gcm -3 . 
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mass of a few x 0.001 M (see Fig. [7J). At high mag- 
netic fields, the distribution is somewhat altered, with no 
material at entropies larger than ~15 &b baryon -1 , and 
ejecta electron fractions extending down to ~0.1. More 
importantly, the total ejecta mass is a factor of ~40 larger 
than in the Bll model and renders such distributions of 
interest. At the end of the B12 simulation, we find a 
total ejecta mass of 0.15 M , with O.lM of the mate- 
rial with IVs between 0.1 and 0.2. and an entropy of 
S ~10/cb baryon -1 . Following the argument of Wheeler 
et al. (1998), in reference to Hoffman et al. (1997), this 
indicates the production of nuclei with an atomic mass of 
up to 190 and motivates further investigation concerning 
r-process nucleosynthesis under such conditions. 

Woosley & Baron (1992) and Fryer et al. (1999) used 
the nucleosynthetic yields of AICs to estimate their oc- 
currence rate. Based on constraints set by a number of 
isotopic abundances (e.g., 62 Ni, 88 Sr; see Woosley & Hoff- 
man 1992), Hartmann et al. (1985) suggested that less 
than 10 -5 Mq of material with an electron fraction below 
0.4 can be ejected per supernova event. By the end of the 
Bll (B12) simulation, a total of -0.004 M (-0.15 M ) 
of such low-F e material has been ejected (note that this 
includes material that has left the grid, the maximum 
computational radius being only 4000 km), which, to- 
gether with the isotopic constraints, sets the rate to a 
maximum of — 5xl0 -5 yr -1 for Bll-type events, and of 
— 1.3xl0 -6 yr -1 for B12-type events. 

Following Woosley & Baron (1992), an alternate con- 
straint can be imposed by using the solar ratio of the 
oxygen to the 88 Sr abundance, assuming oxygen is con- 
tributed primarily by core-collapse explosions of massive- 
star progenitors, and 88 Sr primarily from AICs. For 
the Bll simulation, material with a Y e between 0.4 and 
0.45 is ejected with a very small cumulative mass of 
— 5xl0 -5 M at 800 ms after bounce, corresponding to 
a rate constraint of 0.02 yr _1 . In the B12 simulation, 
more mass is expelled in total, but no ejecta have an 
electron fraction above 0.4 after —200 ms after bounce, 
and we indeed find a total ejected mass with an electron 
fraction between 0.4 and 0.45 that is comparable to that 
of the Bll simulation, and a rate of — 0.02 yr -1 . Hence, 
the second constraint set by the oxygen and 88 Sr solar 
abundances is not really relevant for the AIC we simu- 
late, since we predict nearly all of the nuclei ejected have 
an electron fraction below 0.4. Overall, if the neutrino- or 
magnetically-driven winds we witness in the Bll and B12 
simulations persist over a number of seconds, the corre- 
sponding galactic AIC rate, as limited by nucleosynthetic 
yields, should be at most 9 10 -6 yr -1 , thus from a factor 
of a few to up to two orders of magnitude lower than that 
obtained in the simulations of Woosley & Baron (1992) 
or Fryer et al. (1999). Our B12 simulation suggests 
that if the progenitor core of such AICs rotates fast and 
significant magnetic field amplification occurs after PNS 
formation, the event rate of the AIC of white dwarfs is 
doomed to be very small, of the order of 10 -6 yr -1 , less 
than four orders of magnitude lower than the event rate 
of massive star explosions. 

9 It is unclear how fast the mass loss rate in the B12 model 
is going to decrease over the next few seconds, as it must, and, 
consequently, the corresponding cumulative mass loss cannot yet 
be accurately determined. 



We observe a strong pole-to-equator variation of all 
neutrino fluxes in both the Bll and the B12 models. This 
property was discussed in D06 and is linked to the fast 
rotation and oblateness of the PNS (see also Walder et al. 
2005). Here, we find that this anisotropy persists with 
the introduction of magnetic fields. Moreover, although 
the anti-electron neutrino and the V M " neutrino fluxes 
remain comparable in both the Bll and B12 models, 
we find that the electron neutrino flux is systematically 
lower in the B12 model (for a given direction). 

4.3. Are AICs j -ray-burst /X-ray- flash progenitors? 

The AIC of a white dwarf has been proposed repeat- 
edly as a potential GRB progenitor (Usov 1992,1994; 
Dar et al. 1992; Yi & Blackman 1997,1998; Dai & Lu 
1998). The attractiveness of such collapsed white dwarfs 
stems partly from 1) the absence of a progenitor enve- 
lope, preventing the choking of possibly relativistic ejecta 
launched at the surface of the newly- formed neutron star, 
and 2) an inherently large angular momentum that may 
be insufficient in massive star progenitors to meet the 
conditions necessary for the collapsar model (Woosley 
1993; MacFadyen & Woosley 1999). 

Usov (1992, 1994) invoked the generation of a rel- 
ativistic electron-positron plasma from such a fast- 
rotating highly-magnetized neutron star, powered over 
a timescale of seconds, leading to the emission of X-rays 
and 7-rays at the wind photosphere. A similar model 
was proposed by Yi & Blackman (1997,1998) and Dai 
& Lu (1998), in which a magnetized millisecond pulsar 
powers a fireball by magnetic dipole radiation. Alterna- 
tively, Dar et al. (1992) suggested that electron-positron 
pair production from neutrino-antineutrino annihilation 
outside the PNS would trigger a GRB, provided the neu- 
tron star environment is not polluted with more than 
3xlO _4 M of baryons. 

In the context of standard core-collapse supernova 
explosions, neutrino-antineutrino annihilation does not 
lead to appreciable net energy gain for the wind or ejecta 
because the energy tends to be deposited deep in the 
PNS, or at best in the cooling region in the vicinity 
of the neutrinosphere, thereby contributing negligibly to 
the supernova explosion (Cooperstein et al. 1987). Ad- 
ditionally, the annihilation cross section is strongly en- 
hanced for large-angle collisions, favoring non-spherical 
configurations like disks or tori (Janka 1991; Ruffert & 
Janka 1997). The disk-like morphology of the neutron 
stars formed in the AIC of white dwarfs, and the ad- 
ditional pinching of isodensity contours in the polar re- 
gions, seems to offer an interesting prospect. Following 
the method of Ruffert & Janka (1997) and using the final 
configuration of the B12 model we estimate an v e v e an- 
nihilation power of 0.01 B s _1 , deposited primarily along 
the polar direction, but within 20-50 km of the PNS 
center. This is a rough estimate, but it suggests that 
if our B11/B12 models are sensible representations of 
the AIC of white dwarfs, then little can be gained from 
neutrino-antineutrino annihilation, in particular in com- 
parison with the power injected by magnetic means. This 
is the case at a few hundred milliseconds after bounce, 
but it may change at later times of many seconds. 

In the general context of core-collapse super novae, Tan 
et al. (2001) investigated the acceleration to relativis- 
tic speeds of the surface layers of the progenitor as the 
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shock, generated at bounce, breaks out. The AIC of 
white dwarfs presents an interesting context for such an 
effect, given the absence of a progenitor envelope. But 
the explosion energies and the eject a masses reported 
in the spherical simulations of the AIC of white dwarfs 
(Woosley & Baron 1992; Fryer et al. 1999) were low and 
Tan et al. confirmed the conclusion of those authors that 
the AIC of white dwarfs cannot lead to any sizable 7-ray 
emission or relativistic ejecta. In D06, and reproduced 
here as well for the Bll model, the explosion is strongly 
underenergetic for weak magnetic fields (reaching only 
~0.1B after about 500 ms, with a power that has stabi- 
lized to ~0.1Bs -1 and that may be sustained for a few 
seconds). Most of the power is in the form of a neutrino- 
driven wind, with an associated mass loss rate of a few 
x 0.001 M s _1 . The original blast that breaks out of 
the polar caps of the progenitor white dwarf would rep- 
resent the best option in the context of the Tan et al. 
study, but the shock, debilitated by neutrino losses and 
the photodissociation of nuclei, stalls only mildly. It is 
not re-energized by neutrino energy deposition. Rather, 
it breaks out of the white dwarf by virtue of the abrupt 
density gradient. The resulting early blast amounts to a 
mere 10 -4 B, and despite the low baryon loading of the 
ejecta (~1O _4 M ), they reach velocities of only a tenth 
of the speed of light 10 . In model B12, the explosion en- 
ergy goes above 1 B, but the wind mass loss rate is com- 
mensurate with this energy gain, and is on the order of 
0.1 M s _1 . As a result, the production of a GRB via the 
Tan et al. scenario is again not viable. The asphericity of 
the explosion does not dramatically alter this conclusion, 
since, while the explosion energy is deposited in a smaller 
volume, the mass flux is also confined to a smaller open- 
ing angle and leads to higher effective baryon loading. 
The spherically- symmetric, adiabatic, steady-state, LTE 
model of a super-Eddington wind by Paczynski (1990) 
supports the same conclusion. Given a rate of energy 
deposition of a 0.1-1 B s _1 and a mass loss rate between 
a few x 0.001 and -0.7 M s" 1 for the Bll and B12 
models, respectively, the temperature at the photosphere 
would be vanishingly small in both cases, with a photon 
luminosity of at most 10 40 ergs _1 . The baryon loading 
of the ejecta that occurs in the B12 model, associated 
with the magnetic-driving of the wind, nulls the benefit 
of the additional explosion power. The low photospheric 
temperature obtained here suggests that, since the ejecta 
become optically thin at such a late time, no radiation 
would be visible in the optical. The negligible amount 
of ejecta material with a Y e of 0.5, which corresponds to 
radioactive 56 Ni, means that there is no energy source to 
counteract the strong adiabatic cooling of the neutrino- 
or magnetically-driven wind. It seems that neutrinos and 
gravitational waves would be the only significant signa- 
tures of the explosion. 

The AIC of a white dwarf might give birth to a mag- 
net ar, following magnetic-field compression and ampli- 
fication in the fast rotating neutron star. Yi & Black- 
man (1997,1998), and Dai & Lu (1998) suggested that 
such magnetars could, through extraction of rotational 

10 Note here that the low-density ambient medium surrounding 
the progenitor white dwarf may be slightly inhibiting further ac- 
celeration, but Fig. 11 of Tan et al. 2001 suggests that this shock 
breakout is in any case symptomatically underenergetic. 



energy, power Poynting-flux dominated ejecta and lead 
to GRBs at cosmological distances. Indeed, during the 
cooling phase of such fast-spinning highly-magnetized 
neutron stars, the millisecond proto-magnetar wind may 
transition from non-relativistic and thermally-driven to 
magneto-centrifugally-driven, and finally to relativistic 
and Poynting-flux dominated (Thompson et al. 2004; 
Bucciantini et al. 2006; Metzger et al. 2006; Thomp- 
son 2007). This scenario has in fact the potential to 
make the AIC of white dwarfs detectable in the 7-ray 
range. Furthermore, based on the local Universe dis- 
tribution of short-hard GRBs among all galaxy types, 
rather than exclusively with young stellar populations, 
Chapman et al. (2006) have proposed the AIC of white 
dwarfs as an alternative to massive star progenitors for 
the formation of magnetars, whose soft 7-ray emission is 
associated with giant-flaring events. We concur in prin- 
ciple with this proposition, although the AIC rate of at 
most ~10 -6 yr _1 for B12-type events that we infer in §4.21 
suggests that magnetars originating from the AIC of a 
fast-rotating white dwarf should be significantly under- 
represented in the magnetar population. 

An AIC analogue to the collapsar model (Woosley 
1993; MacFadyen & Woosley 1999) is an obvious alterna- 
tive, with both the transition of the PNS to a black hole 
and energy deposition in the excavated polar regions at a 
high rate of a few x 0.1 B s _1 . In the context of the AIC 
of white dwarfs, the potential for formation of a black 
hole seems uncertain, as it might be difficult in nature 
to ultimately meet the basic mass requirement of a 2- 
2.5 M (baryonic) neutron star for black- hole formation 
(Baym et al. 1971; Arnett & Bowers 1977). In D06, we 
studied two models, with 1.46 M and 1.92 M masses. 
The PNS mass was larger in the lower mass progenitor 
because in the 1.92 M model much of the mass ended in 
a massive disk. Higher-mass AIC progenitors would need 
even more angular momentum to achieve a stable equi- 
librium configuration, which would likely result in the 
formation of even higher-mass disks on extended quasi- 
Keplerian orbits. At the same time, magnetic stresses at 
the surface of the neutron star lead to such high mass 
loss rates that the PNS mass paradoxically decreases 
with time (the PNS mass in the B12 model decreases 
to ~1.3M at ~500ms after bounce). It, thus, seems 
that accretion, leading to black-hole formation, would 
have to operate over a very extended time, and it is not 
clear that given the explosion power that can be obtained 
(witness the results for the B12 model) whether accre- 
tion will ever win against explosion. A further obstacle 
to high accretion rates in the AICs is the extended con- 
figuration of the mass distribution outside the neutron 
star, with progenitor lobes located at ~1000km, rather 
than around ^100 km in the ideal setup for the collapsar 
model (Woosley 1993; MacFadyen & Woosley 1999). The 
outcome here is closely related to the original mass distri- 
bution of the white dwarf at the time of collapse. Current 
SPH simulations of coalescing white dwarfs, one forma- 
tion channel of AICs, suggest that a significant fraction 
of that material goes into a Keplerian disk, loosely-bound 
and, thus, not directly involved in the collapse/bounce 
of the core, nor with the subsequent explosion. This 
material may represent a substantial amount of mass, 
perhaps up to 1 M , but is located at least at 1000 km 
away from the neutron star or would-be black hole, so 
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that the actual accretion rate is very small. In the con- 
figuration of the 1.92 M model described here, we have 
0.5 Mq available in the disk, with an associated viscous 
accretion timescale of ~100s (using an a parameter of 
0.1; Shakura & Sunyaev 1973), and hence, a feeble accre- 
tion rate of 5xlO _3 M0 s _1 . This is a factor of a hun- 
dred lower than the minimum value quoted in the con- 
text of the collapsar model. This extended configuration 
for AICs is in fact reminiscent of the fallback scenario of 
Rosswog (2007) for double neutron star mergers, wherein 
material is placed on eccentric orbits during the merger 
event. Note that, here again, accretion onto a black hole 
is a prerequisite for the powering of a relativistic jet, and, 
ultimately, 7-ray radiation. In the context of AIC, such 
accretion may simply extend the PNS-wind phase over a 
longer period, and not lead to the formation of a black 
hole. 

4.4. Ultimate fate of the AIC of white dwarfs 

The AIC of white dwarfs seems very unfavorable as 
a progenitor of 7-ray, or even optical bursts, and direct 
detection of such events might instead come from neu- 
trino or gravitational wave detectors (see ^3] and §10 of 
D06). Indirectly, we may infer their existence through 
the r-process nuclei ejected during the explosion ( M4.2p . 

Another avenue is to search for them in the neutron 
star population. Fallback following SN explosions of mas- 
sive progenitors is expected to bury the PNS and weaken 
its surface magnetic field 11 , as in magnetic neutron stars 
accreting from a companion star (see, e.g., Geppert & 
Urpin 1994; Payne & Melatos 2004,2007). This argument 
would suggest that magnetic neutron stars may originate 
from the AIC of white dwarfs, since these are exempt 
from such fallback. Observations of magnetic neutron 
stars in binary systems have indeed been interpreted in 
this way (see, e.g., Taam & van den Heuvel 1986, Ergma 
1993; van den Heuvel & Bitzaraki 1995; van Paradijs et 
al. 1997), although our simulations suggest that some 
"fallback", albeit modest, could occur from the progen- 
itor lobes of the white dwarf and bury the large surface 
fields of ~10 15 G that we see in the B12 model a few hun- 
dred milliseconds after bounce. Hence, in this context, 
whether additional material from the progenitor lobes 
of the white dwarf accretes onto the highly-magnetized 
PNS may determine whether it will retain its magnetar- 
like fields or transition to a weaker-field pulsar. Despite 
the extraction of angular momentum by magnetic fields, 
the PNS remains a millisecond-period rotator. Further 
accretion from the Keplerian disk would mitigate the loss 
of angular momentum, and, thus, AICs of white dwarfs 
should consistently lead to rapidly-rotating pulsars af- 
ter this original phase. On longer timescales of hours or 
days, however, magnetic torquing may be more in evi- 
dence and may spin down the PNS as it transitions to a 
magnetar-like object. In this case, the remnant could be 
a highly-magnetized object, with a much reduced rota- 
tion rate, a correlation that has been found in the pulsar 
population (Vranevsevic et al. 2007). Interestingly, if 
such accretion is not significant, this material could be 

11 In this context, the mass accreted after a delay starts with a 
very weak initial field. Little compression ensues and the spin up 
is moderate, so that, despite the presence of velocity shear in the 
region of accretion, we expect only modest field amplification. 



found in a debris disk, as recently observed around an 
isolated neutron star (Wang et al. 2006). 

Finally, we reported in D06 that in the simulation for 
the 1.46 M AIC model of Yoon & Langer (2005), car- 
ried over 180° (i.e. including both hemispheres), perfect 
north-south symmetry was preserved. We have used this 
observation to carry out all subsequent investigations in 
a 90° quadrant, thus in just one hemisphere. We expect 
that the AIC of white dwarfs leads to essentially perfect 
north-south symmetry, a property that applies more gen- 
erally to the explosion following the collapse of ONeMg 
cores, universally characterized by a steep density profile 
above the core. Consequently, these explosions should 
involve no significant kick to the PNS. This supports 
the notion that the AIC of white dwarfs can, at least 
theoretically, contribute to the neutron star population 
of globular clusters (see, e.g., Bailyn & Grindlay 1990; 
Ivanova et al. 2006). But given the very low inferred 
rates of AICs (in particular for the B12-like events), the 
ONeMg cores of low-mass massive stars are more likely 
candidates. Similarly, this also supports the association 
of collapsed ONeMg-cores with the populations of low- 
velocity, low-mass neutron stars found in low-eccentricity 
double neutron star systems (see, e.g., van den Heuvel 
2007). 

5. CONCLUSION 

We have presented multi-group, flux-limited diffusion, 
magnetohydrodynamics simulations of the Accretion In- 
duced Collapse (AIC) of white dwarfs, focusing on the 
1.92 M model of Yoon & Langer (2005). The main in- 
centives were to gauge the relevance and estimate the 
magnitude of MHD effects in such collapse events, and 
thereby to address a potential limitation of the simula- 
tions in D06. The key ansatz in the present work is that, 
although the magnetic field distribution and magnitude 
in the progenitor white dwarf are unknown, fast and dif- 
ferential rotation in the collapsed object is expected to 
lead to an exponential amplification of any initially weak 
poloidal field through the magnetorotational instability, 
reaching large values at saturation. In this context, the 
magnitude of the post-bounce MHD effects we report 
from our simulations will thus be determined, in a more 
physical sense, by the angular momentum budget rather 
than by the initial field in the pre-collapsed white dwarf. 
Provided angular momentum is abundant, we find that 
MHD effects lead to large quantitative changes when 
compared with an equivalent situation ignoring magnetic 
fields. Said differently and in a broader sense, we expect 
significant MHD effects in core-collapse events with fast 
rotation, effects that cannot be ignored without compro- 
mising the conclusions 12 . 

Perhaps the most striking property of the AIC of a 
1.92 M white dwarf model with large initial magnetic 
fields is that the explosion energy at just a few hundred 
milliseconds after bounce is no longer an anemic 0.01 B, 
but is on the order of one Bethe, with an ejecta mass 
raised from a few x 0.001 M to ~O.1M . Significant 

12 Non-axisymmetric instabilities of such fast-rotating white 
dwarfs may alter the angular momentum distribution after bounce 
(Ott et al. 2006), with the potential to modify magnetic-field am- 
plification rate as well as the accretion rate we observe in our sim- 
ulations. Numerical simulations of the kind presented here but in 
3D are needed to quantify these effects. 
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rotational energy is extracted from the core with the in- 
clusion of magnetic fields, directly seen in the decrease in 
rotational energy and the spin down of the PNS. This en- 
ergy is transformed into magnetic energy that powers the 
explosion. For higher B-fields, the density distribution 
at the surface of the neutron star is more spherical, the 
10 10 gcm -3 and 10 11 gcm -3 density contours stretching 
further out along the pole, aided by magnetic stresses. 
We find that mass loading of the eject a is so severe that 
no relativistic high-energy electromagnetic transient can 
be expected during the PNS phase. In the past, the AIC 
of white dwarfs has been invoked as a potential GRB 
progenitor, but this was usually done in ignorance of the 
neutrino-driven wind that prevails during the first tens 
of seconds after the formation of the PNS. MHD effects 
seem to make this situation even worse, not better, by 
loading magnetically-driven ejecta with 100 times more 
baryons. 

In the single-degenerate formation scenario, magnetic 
torques operating over evolutionary timescales might 
force the white dwarf into solid-body (and, perhaps, crit- 
ical) rotation, and prevent its growth much beyond the 
canonical Chandrasekhar mass, thus preventing a tran- 
sition to a black hole after PNS formation. However, 
in the double-degenerate scenario, differential and fast 
rotation results after the merging event and a super- 
Chandrasekhar white dwarf may form. Ultimately, tran- 
sition to a black hole may occur if the total PNS mass 
exceeds the general-relativistic limit for gravitational sta- 
bility. The transition to a black hole, by shutting off the 
wind from the PNS, may then lead to a 7-ray burst or an 
X-ray flash, but the 2-2.5 M mass limit for black hole 
formation restricts the number of adequate AIC progen- 
itors to a very small number, difficult to estimate, but 
perhaps to one in a hundred AICs. The coalescence of 
two white dwarfs 13 seems the more favorable configu- 
ration, but no more than a few x 0.1]Vi©of material 
placed at large distances (few 1000 km) would be avail- 
able for accretion, and this over a viscous timescale of 
~100s and, thus, with a feeble accretion rate of at most 
5xlO _3 M s _1 . In the unlikely event that all these 
necessary conditions were met, the relativistic beaming 
(with opening angle 0- } ) of the radiating ejecta would 
decrease the probability of 7-ray detection by a factor 
0.5(1 — cos#j), i.e., 2-3 orders of magnitude for a 5°-15° 
opening- angle 14 . Together with the constraints from 
nucleosynthetic yields of B12-like events, the 1 in a 100 
AICs that may reach the black-hole mass limit, and the 
inherent beaming of the resulting 7-ray radiation, we es- 
timate the detection rate of 7-ray bursts associated with 
AICs to be on the order of 10 -10 yr _1 , significantly lower 
than the average rate per galaxy of 4 x 10 -7 yr _1 inferred 
from GRBs detected with BATSE (Zhang & Meszaros 
2004). 

Moreover, the absence of nickel in the ejecta precludes 

13 It even seems that these two white dwarfs would have to 
have comparable masses. Otherwise, the secondary would likely 
get disrupted and with its matter dispersed in a disk around the 
primary, rather than merging with the primary (Guerrero et al. 
2004). 

14 This relativistic beaming of 7-ray radiation is inferred from 



any significant optical signature for the event, which 
should thus remain obscure at all wavelengths. If the AIC 
occurs in our Galaxy, possible signatures should instead 
be from the neutrino burst that accompanies the forma- 
tion and the cooling of the PNS, and gravitational waves 
from the very aspherical bounce and the anisotropic neu- 
trino emission (Ott 2007). 

The AIC of white dwarfs may, however, have a sig- 
nificant impact on the interstellar medium composition 
in r-process nuclei. We find that the AIC of rapidly- 
rotating white dwarfs leads to explosions that eject co- 
pious amounts of neutron-rich material, with an electron 
fraction reaching down to ~0. 1-0.2, depending on the 
magnitude of the neutrino/anti- neutrino absorption by 
the ejecta. If not magnetically-aided, the neutrino- driven 
wind leads, however, to a low ejecta mass of a few x 
0.001 M , which, given the low occurrence rates for AIC, 
makes these yields somewhat irrelevant. If magnetically- 
aided, as in the simulation B12 presented here, about 
0.1 M of material is ejected with an electron fraction be- 
tween 0.1 and 0.2, making the AIC of white dwarfs a pos- 
sible r-process site, as promoted in an analogous context 
by Wheeler et al. (1998). More importantly, constraints 
placed by solar isotopic abundances require that such 
pollution of the interstellar medium by magnetically- 
driven winds of AICs occurs with an event rate of at 
most ~10 _6 yr _1 . If we believe the context adopted for 
this work, this suggests that the circumstances leading 
to such fast rotating white dwarfs are rare. 

An important ramification of this work is that, free of 
significant fallback, the AIC of white dwarfs offers a pos- 
sible scenario for the formation of highly magnetized neu- 
tron stars, potentially with magnetar-like surface fields of 
10 15 G, with magnetically- tor qued cores that spin down 
over many seconds. Here again, though, the low event 
rate of such AICs should result in their contributing only 
modestly to the magnetar population. 
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the achromatic steepening breaks of some GRB optical afterglow 
lightcurves (Kulkarni et al. 1999; Harrison et al. 1999). 
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